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Abstract 

The present paper is concerned with the development of a micromechanical model of the harden- 
ing, rate-sensitivity and thermal softening of bcc crystals. In formulating the model we specifically 
consider the following unit processes: double-kink formation and thermally activated motion of kinks; 
the close-range interactions between primary and forest dislocations, leading to the formation of jogs; 
the percolation motion of dislocations through a random array of forest dislocations introducing short- 
range obstacles of different strengths; dislocation multiplication due to breeding by double cross-slip; 
and dislocation pair annihilation. The model is found to capture salient features of the behavior of Ta 
crystals such as: the dependence of the initial yield point on temperature and strain rate; the presence 
of a marked stage I of easy glide, specially at low temperatures and high strain rates; the sharp onset of 
stage II hardening and its tendency to shift towards lower strains, and eventually disappear, as the tem- 
perature increases or the strain rate decreases; the parabolic stage II hardening at low strain rates or high 
temperatures; the stage II softening at high strain rates or low temperatures; the trend towards saturation 
at high strains; the temperature and strain-rate dependence of the saturation stress; and the orientation 
dependence of the hardening rate. 
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1 Introduction 



The present paper is concerned with the development of a micromechanical model of the hardening, rate- 
sensitivity and thermal softening of bcc crystals. We place primary emphasis on the derivation of closed- 
form analytical expressions describing the macroscopic behavior of the crystals amenable to implementation 
as constitutive relations within a standard finite-element code. In developing the model, we follow the well- 
established paradigm of micromechanical modeling, consisting of: the identification of the dominant or rate- 
limiting 'unit' processes operating at the microscale; the identification of the macroscopic forces driving the 
unit processes; the analysis of the response of the unit processes to the macroscopic driving forces; and the 
determination of the average or macroscopic effect of the combined operation of all the micromechanical 
unit processes (see, e. g., [Jl], ||, ||, ||, ^ for recent overviews and discussions pertaining to micromechanics 
and multiscale modeling of materials). 

In formulating the present model we specifically consider the following unit processes: double-kink 
formation and thermally activated motion of kinks; the close-range interactions between primary and for- 
est dislocations, and the subsequent formation of jogs; the percolation motion of dislocations through a 
random array of forest dislocations introducing short-range obstacles of different strengths; dislocation mul- 
tiplication due to breeding by double cross-slip; and dislocation pair annihilation. We believe that this 
forms an 'irreducible' set of unit processes, in that each of these processes accounts for-and is needed for 
matching-salient and clearly recognizable features of the experimental record. For instance, consideration 
of dislocation multiplication and annihilation leads to a predicted stage I-II transition strain which decreases 
with temperature and increases with strain rate, as observed experimentally. 

We bring an assortment of analysis tools to bear on each of the unit processes under consideration. 
As already mentioned, the choice of tools is to a large extent conditioned by our desire to derive closed- 
form analytical expressions for all constitutive relations. The motion of dislocations through an otherwise 
defect-free lattice is assumed to be thermally activated and controlled by the nucleation of kink pairs. A 
analysis of this process based on transition-state theory yields the effective Peierls stress of the lattice as 
a function of temperature and strain rate. The short-range interactions between a primary and secondary 
dislocation are assumed to result in the acquisition by both dislocations of a jog. The process by which the 
primary dislocation unpins from a forest obstacle is assumed to be thermally activated and the activation 
energy barrier is identified with the jog-formation energy. These assumptions result in temperature and 
rate-dependent obstacle strengths which are in good qualitative and quantitative agreement with experiment. 

The percolation motion of dislocations through random arrays of point obstacles is studied in detail 
using statistical methods. The analysis presented here generalizes the analysis of Cuitino and Ortiz [Q] 
so as to account for multi-species distributions of finite-strength obstacles. Finally, we model dislocation 
multiplication as the result of two competing effects: the proliferation of dynamic sources by double cross- 
slip; and dislocation pair annihilation. We model cross-slip as a thermally activated process and develop 
an elastic model, similar to that of Huang et al. [§], which estimates the probability that the trajectories 
of two parallel screw dislocations collide, resulting in the annihilation of the pair. The model predicts that 
multiplication and annihilation balance out at sufficiently large strains, leading to saturation. 

We validate the micromechanical model by recourse to detailed comparisons with the uniaxial tension 
tests on Ta single crystals of Mitchell and Spitzig The model is found to capture salient features of 
the behavior of Ta crystals such as: the dependence of the initial yield point on temperature and strain rate; 
the presence of a marked stage I of easy glide, specially at low temperature and high strain rates; the sharp 
onset of stage II hardening and its tendency to shift towards lower strains, and eventually disappear, as the 
temperature increases or the strain rate decreases; the parabolic stage II hardening at low strain rates or 
high temperatures; the stage II softening at high strain rates or low temperatures; the trend towards satura- 
tion at high strains; the temperature and strain-rate dependence of the saturation stress; and the orientation 
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dependence of the hardening rate. 



2 General Framework 

Our objective is to formulate a model of the hardening of bcc crystals which is well-suited to finite- 
deformation large-scale finite-element calculations and, therefore, we couch the model within a nonlinear 
kinematics framework. To this end, we adopt a conventional multiplicative elastic-plastic kinematics of the 
form 0, 111], 0, |B[] 



F = F e F p (1) 

where F is the deformation gradient, F e is the elastic lattice distortion and rotation, and F p is the plas- 
tic deformation, which represents the net effect of crystallographic slip and leaves the lattice undistorted 
and unrotated. In addition, the plastic deformation is volume preserving. In materials such as metals, the 
elastic response is ostensibly independent of the internal processes and the free energy density decomposes 
additively as 

A = W e (F e ,T) + W P (T,F P ,Q) (2) 

where T is the absolute temperature and Q is some suitable set of internal variables. The function W e 
determines the elastic response of the metal, e. g., upon unloading, whereas the function W p describes the 
hardening of the crystal. Physically, W p measures the energy stored in the crystal due to the plastic working 
of the material. 

The examples of validation reported subsequently probe the hardening behavior of the material due 
to crystallographic slip over a range of moderate temperatures and strain rates. We therefore eschew a 
number of issues which become important at high pressures and elevated temperatures, but which play a 
limited role otherwise, including the volumetric equation of state (EoS), the pressure dependence of yield 
and elastic moduli, and others. First-principles calculations of the EoS and elastic moduli of bcc metals up 
to high pressures and temperatures may be found in [16, [n|, 18, 19, 21]. These results may be used as 



a basis for extending the present theory to conditions such as arise in shocked materials. These extensions 
notwithstanding, a simple form of the elastic energy density appropriate for present purposes is 

W e (F e ,T) = X - Cl]kl (T){e% - a l3 T)(e e kl - a kl T) (3) 
where c is the tensor of elastic moduli, 

1. 



e e = logVC e = -log(F ei F e ) (4) 

is the logarithmic elastic strain, and a is the thermal expansion tensor. The elastic properties of a cubic 
crystal are fully described by the three Voigt constants en, c\2 and C44 and a scalar thermal expansion 
coefficient a. For low temperatures, the temperature dependence of the moduli may be taken to be linear in 
T to a first approximation, leading to the relation: 

Ci jk i(T) & c1 jkl - T cj jhl (5) 



This linear dependence of the elastic moduli is observed experimentally up to moderate temperatures y22|], 
and can be justified within the framework of statistical mechanics 
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Table 1 : Slip systems of bcc crystals 



Plastic deformations in single crystals are crystallographic in nature. Following Rice [15], we adopt a 
flow rule of the form 



F P = ^ j a s a <g> m a ^j F p (6) 

where 7" is the slip strain, and s a , and m a are orthogonal unit vectors defining the slip direction and slip- 
plane normal corresponding to slip system a. The collection 7 of slip strains may be regarded as a subset 
of the internal variable set Q. A zero value of a slip rate 7° signifies that the corresponding slip system 
a is inactive. The flow rule @ allows for multiple slip, i. e., for simultaneous activity on more than one 
system over a region of the crystal. The vectors {s a , m a } remain constant throughout the deformation and 
are determined by crystallography. For bcc crystals, we consider the 24 slip systems listed in Table [I] in the 
Schmid and Boas nomenclature. 

Slip on the {112} systems is known to be asymmetric at low temperatures: slip is easier when the 



applied stress is such that dislocations move in the twinning direction (e. g., [24, 25, 26]). However, over 
the temperature range of interest here the experimental evidence [26, |7|] suggests that the extent of the 



twinning/anti-twinning asymmetry is relatively small, e. g., of the order of 20 MPa for Mo at 150 K [26], 
and, for simplicity, we will neglect it to a first approximation. 

In the examples presented in Section [7|, the constitutive equations are integrated in time using the vari- 
ational update of Ortiz and Stainier [E|]. The variational formulation of the rate problem proposed by 



Ortiz and Stainier [ ]28[ ] additionally furnishes a convenient avenue for the superposition of the various unit 
mechanisms analyzed subsequently. Our working assumption is that the dissipation rates arising from these 
mechanisms are additive. This assumption in turn implies that the critical resolved shear stress r" for the 
operation of system a may be computed as the sum of unit-process contributions, i. e., 

r c Q = r" (process) (7) 



processes 



where the sum extends to all unit processes. In the present model, these contributions are the Peierls stress 
of the lattice and the forest-obstacle resistance. Other dissipation mechanisms, such as phonon drag, may in 
principle be superposed likewise. 
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3 Dislocation mobility 



In this section we begin by considering the thermally activated motion of dislocations within an obstacle- 
free slip plane. Under these conditions, the motion of the dislocations is driven by an applied resolved 
shear stress r and is hindered by the lattice resistance, which is weak enough that it may be overcome by 
thermal activation. The lattice resistance is presumed to be well-described by a Peierls energy function, 
which assigns an energy per unit length to dislocation segments as a function of their position on the slip 
plane. 

In bcc crystals, the core of screw dislocation segments relaxes into low-energy non-planar configurations 
[ p^ , 30, 31, 32|, [?3|, ||, 34, 35]. This introduces deep valleys into the Peierls energy function aligned with the 



Burgers vector directions and possessing the periodicity of the lattice. At low temperatures, the dislocations 
tend to adopt low-energy configurations and, consequently, the dislocation population predominantly con- 
sists of long screw segments. In order to move a screw segment normal to itself, the dislocation core must 
first be constricted, which requires a substantial supply of energy. Thus, the energy barrier for the motion of 
screw segments, and the attendant Peierls stress, may be expected to be large, and the energy barrier for the 



motion of edge segments to be comparatively smaller. For instance, Duesbery and Xu [ ]36[ ] have calculated 
the Peierls stress for a rigid screw dislocation in Mo to be 0.022/i, where \i is the (111) shear modulus, 
whereas the corresponding Peierls stress for a rigid edge dislocation is 0.006//, or about one fourth of the 
screw value. This suggests that the rate-limiting mechanism for dislocation motion is the thermally activated 



motion of kinks along screw segments ([]37|, |38| , p9|]). 

Consider now a screw dislocation segment, possibly pinned at both ends by an obstacle pair, subjected 
to a resolved shear stress r > 0. At sufficiently high temperatures a double-kink may be nucleated with the 
assistance of thermal activation (e. g., pOj [flj , ^J, and the subsequent motion of the kinks causes the screw 
segment to effectively move forward, Fig. 1. Neglecting entropic effects and invoking Peach- Koehler's 
formula, the forward and backward activation enthalpies AG + and AG~, respectively, for the formation of 
a double kink may be approximated as 

AG ± « £ kink =F brL kink l P (8) 

where E is the energy of formation of a kink-pair, L is the length of an incipient double kink, and 
lp is the distance between two consecutive Peierls valleys. For bcc crystals, lp = y/2 /3a if the slip plane 
is {110}, lp = y/2a, if the slip plane is {112}, and lp = \f&j?>a if the slip plane is {123}, where a is the 
cubic lattice size p6[]. 

The formation energy E kmk cannot be reliably estimated from elasticity since is is composed mostly 
of core region. It can, however, be accurately computed by recourse to atomistic models. For instance, for 



Mo at zero stress Xu and Moriarty [41] have found formation energies E kmk of the order of 1 eV for kinks 
separated by a distance greater than L kmk = 156. The core structure, gamma surfaces, Peierls stress, and 
kink-pair formation energies associated with the motion of a/2 (111) screw dislocations in Ta and Mo []|] 
have also been calculated by Moriarty et al. 

We may expect a proliferation of kinks at a stress of the order of 

j^kink 

r ° = bL^Hp (9) 



for which AG = 0. For Mo, Xu and Moriarty [41] have computed tq to be of the order of a few GPa. The 



activation free enthalpy for double kink nucleation may thus be rewritten as 



ACVMri t: KwK [ lT^r ) ■ <>£ r < m (10) 
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A slightly more elaborate empirical formula that is widely used to fit activation energies was proposed by 



Kocks et al. [ |42[ ] and the corresponding parameters have been determined for Ta by Tang al. U43Q, but this 
enhancement will not be pursued here. 

Transition state theory predicts that a dislocation segment effects vd&~^ G and vp,eT^ G jumps per unit 
time in the positive and negative directions, respectively. Here, j3 = 1/kpT, kp is Boltzmann's constant, 
T is the absolute temperature, and vp> is the attempt frequency which may be identified with the Debye 
frequency to a first approximation. Since the length of the individual jumps is lp, the mean velocity of the 
dislocations then follows as 

v = 2lpv D e-? Ekink sinh (W^) dD 
An application of Orowan's formula then gives the slip strain rate as 

j = % kink e -/^«* ginh (pE^l) (12) 



where 



7 kink = 2bpl P u D (13) 



is a reference strain rate, and p is the dislocation density. In writing (|l2|), we have taken into account 
the possibility of thermally activated jumps in the reverse direction, with forward and backward activation 
enthalpies AG + and AG~, respectively. For slip in the positive direction, 7 > 0, ( |I2| ) may be inverted to 
give: 

^ = 1 asinh ( _J_ e /^ kink ^ ( 14) 



T f3E kmk V7o 

where rp may be regarded as a temperature and rate-dependent effective Peierls stress. 

Figure 2 illustrates the dependence of the effective Peierls stress on temperature and rate of deformation. 
The Peierls stress decreases ostensibly linearly up to a critical temperature T c , beyond which it rapidly 



tends to zero. These trends are in agreement with the experimental observations of Wasserbach [ ]44[ | and 
Lachenmann and Schultz [pq]. The critical temperature T c increases with the strain rate. In particular, 
in this model the effect of increasing (decreasing) the strain rate has an analogous effect to decreasing 



(increasing) the temperature, and vice- versa, as noted by Tang et al. [|46p . In the regime of very high strain- 
rates (7 > 10 5 s~ l ), effects such as electron and phonon drag become important and control the velocity of 
dislocations @B8t]. 



4 Forest hardening 



In the forest-dislocation theory of hardening, the motion of dislocations, which are the agents of plastic 
deformation in crystals, is impeded by secondary -or 'forest'- dislocations crossing the slip plane. As 
the moving and forest dislocations intersect, they form jogs or junctions of varying strengths [ p9j |50| , [|, 
51, 52, 53, 8, 54, which, provided the junction is sufficiently short, may be idealized as point obstacles. 



Moving dislocations are pinned down by the forest dislocations and require a certain elevation of the applied 
resolved shear stress in order to bow out and bypass the pinning obstacles. The net effect of this mechanism 
is macroscopic hardening. Detailed numerical simulations of a dislocation line propagating through forest 



dislocations have been carried out by Foreman and Makin [56, 57], and by Kocks [58], and more recently 
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by Tang et al. [ftqj. Analytical treatments of the model have been given by Kocks [pq], Ortiz and Popov 
[p^l, and Cuitino and Ortiz [Q] for the case of infinitely strong obstacles. A phase-field model of the forest 
hardening mechanism has been proposed by Ortiz [[K]] and by Cuitino et al. [El]]. Here we extend the 
statistical analysis of Cuitino and Ortiz ^ to the case of several species of obstacles of finite strength. 
Because of the random nature of the interactions, the motion of dislocations through a distribution 



of obstacles is best described in statistical terms [59]. We begin by treating the case of infinitely strong 



obstacles. In this case, pairs of obstacles pin down dislocation segments, which require a certain threshold 
resolved shear stress s in order to overcome the obstacle pair. Since the distribution of point obstacles 
within the slip plane is random, it follows that s is itself a random variable. We shall let /^(s, t) denote the 
probability density function of two-point barrier strengths on slip system a at time t. The time dependence 
of /^(s, t) is a consequence of the variation in forest dislocation density. 

In order to determine the precise form of Jq(s, t), we begin by noting that the Peierls energy landscape 



of bcc crystals strongly favors either screw or edge segments [52]. In addition kinks, or points of change 
of direction of the dislocation line, carry a non-negligible amount of energy. It therefore follows that the 
lowest-energy configuration of unstressed dislocation segments spanning an obstacle pair is a step of the 
form shown in Fig. 3. Under these conditions, the bow-out mechanism by which a dislocation segment 
bypasses an obstacle pair may be expected to result in the configuration shown in Fig. 3 (bold line). If 
the edge-segment length is l e , a displacement da e of the dislocation requires a supply of energy equal to 
2U screw da e +bTp ge l e da e in order to overcome the Peierls resistance t p ge and to extend the screw segments. 
The corresponding energy release is brl e da e . Similar contributions result from a displacement da s of the 
screw-segment of length l s . Therefore, for the bow-out of the dislocation to be energetically possible we 
must have 

brl e da e + brl s da s > 2U SCTew da e + bT e p dge l e da e + 2U edge da s + bT s P cmw l s da s (15) 

As already noted, in bcc crystals the core energy [/ scrcw per unit length of the screw segments is smaller than 
the core energy U edgc per unit length of the edge segments (e. g., [|3^]). Conversely, the Peierls resistance 
r jcrew to (jj e glide f scr ew segments is larger than the Peierls resistance rp dge for edge segments (e. g., 



[pq]). Retaining dominant terms only, Eq. (|15j) simplifies to 



r > 8 = rf> crew + (16) 

01 s 



We may further identify the Peierls stress r|, crew with (14), whence it follows that s comprises a thermally- 
activated and rate-sensitive term rp and an athermal and rate-insensitive term 2U edge /bl s . 

We shall assume that the point obstacles are randomly distributed over the slip plane with a mean density 
n a of obstacles per unit area. We shall also assume that the obstacle pairs spanned by dislocation segments 
are nearest-neighbors in the obstacle ensemble. Thus, if r is the distance between the obstacles in an ob- 
stacle pair, then the circle of radius r centered at either obstacle contains no other obstacles. Under these 
conditions, the probability density of r is given by [|58|, 0] 

fo(r,t) = 2irn a rex.p(-irn a r 2 ) (17) 

where the dependence of /o on time stems from the time-dependence of the obstacle density n a . 

In order to deduce the probability density of l s , and by extension of s, we note that, given a point obstacle 
located at the origin, the probability of finding another obstacle in the element of area rdrdO is: 

dFft(r,t) = n a rexp(-TTn a r 2 )drd9 (18) 
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Changing variables to Cartesian coordinates, x = r cos 9,y = r sin 0, ( |18| ) may be recast in the form: 

dPg = n a exp[-vrn a (x 2 + y 2 )]dxdy (19) 
The frequency of obstacle pairs with screw-segment length l s is, therefore: 

fo(l s ,t) = 2 n a exp[-7rn Q (x 2 + l 2 s )}dx = Vn« exp(-Tm a l 2 ) (20) 



The corresponding probability of finding an obstacle pair of strength s follows from (2C) and (|16|) as: 



/oM) 



b{s-T«y exp 



7rn 
l 2 " 



2 [/edge 



S — Tr 



S > Tp 



and the associated distribution function is: 



1 — erf 



edge 



b(s - T«) 



(21) 



(22) 



We note that the Peierls stress Tp depends on the slip system a through its dependence on the slip-strain rate 
7 a . 

It is interesting to note that the probability density Jq(s, t) of obstacle-pair strengths just derived for bec 
crystals differs markedly from those which are obtained for fee crystals by a similar argument [58, 63, 64|,f7|], 
namely: 



27m a U 2 
b 2 s 3 



exp 



n a TT U z 



(23) 



where U is the dislocation core energy per unit length of dislocation. This difference owes to the different 
bow-out configurations for the two crystal classes and the comparatively larger values of the Peierls stress in 



bec crystals. Thus, the Peierls stress of fee crystals is generally quite small and, as in the derivation of (23) 
is often neglected entirely to a first approximation. For bec crystals, the effective Peierls stress Tp decreases 



with increasing temperature and decreasing rate of deformation, Eq. (14), and hence the behavior of bec 
crystals may be expected to be closer to that of fee crystals under those conditions, as noted by Tang et al. 
[p^]. Conversely, the hardening behavior of bec crystals may be expected to differ sharply from that of fee 
crystals at low temperatures and high rates of deformation. 

The function /q*(s, t) just derived provides a complete description of the distribution of the obstacle-pair 
strengths when the point obstacles are of infinite strength and, consequently, impenetrable to the disloca- 
tions. Next we extend the preceding analysis to the case of finite obstacle strengths. Let s al3 be the strength 
of the jogs or junctions formed by dislocations of systems a and j3. Now consider an obstacle pair in which 
the weakest point obstacle corresponds to a forest dislocation of type (3 and, therefore, has strength s Q/3 . The 
probability that the strength of the obstacle pair be s is, therefore, 



1-H(s- s aP ) 



(24) 



where H(s) is the Heaviside function. The probability that an obstacle on system a be of type (3 ^ a is 
n a @ /n a , where n a ^ is the number of obstacles of type (3 per unit area of the slip plane a, and 



n 



*(t) = 5>^(t) 



(25) 
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is the total obstacle density on slip system a. 

Next we note that the probability that the weakest of the two obstacles forming a obstacle pair be of type 

(3 is: 



P 



It is readily verified that 



-.a/3 I „af3 
a(3 _ " 1 " 



- + 2 y 



n 



a/3' 



s af3' >s af3 



n u 



, (3 ^ a, f3 7^ a 



(26) 



y P °p 

(3^a 



.a/3 

T( — 



(3^a 



(27) 



as required. Finally, the probability that an obstacle pair have strength s follows as 

r( S ,t) = y P ^(t)r( S \s^,t) 

(3^a 



(28) 



This probability distribution function jointly accounts for the strength of the obstacle pairs due to line tension 
and to the obstacle strength. Evidently, f a (s,t) is supported in the interval [0,s max ], where s max is the 
maximum obstacle strength. Making use of (Q), the probability density function (28) can be rewritten in 
the form 



ns,t) = y P ^(t)fZ(s,t) 



(3^a 

and the associated distribution function becomes 



1 - H(s - s afj ) 



(29) 



f3+a I 



l-H{s- s 



a/3) 



+ H(s - s' 



a/3\ 



(30) 



We assume that n al3 , the number of obstacles of type (3 per unit ai - ea of the slip plane a, scales with the 
dislocation density p@ according to the relation: 



(3D 



where the coupling constants a ap are regarded here as purely geometrical parameters. A simple geometrical 
argument based on counting intersections of randomly distributed lines with a slip plane gives [ 65 ] : 



n 



a/3 



7T 



1 - (m a ■ mP) 2 pP 



(32) 



where m a and mP are the unit normals to slip planes of types a and and p' 3 is the dislocation line density 
per unit volume in slip system j3. A comparison between ( pi] ) and (|2|) suggests writing: 



,a/3 



7T 



ao-y 1 - (m a ■ mf 3 ) 2 



(33) 



where ao < 1 is an 'efficiency' factor which accounts for the tendency of dislocations to tangle and form 
loops, which in turn tends to lower the number of slip-plane crossings. 
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Let T a (t) now be the resolved shear stress acting on the slip system a at time t. Assume for now 
that r a (t) increases monotonically from zero at t = 0. Evidently, for dislocations to be stable at time t 
they must face barriers of strengths s in excess of T a (t). As r a (t) is increased to r a (t) + f a (t)dt, the 
dislocation segments held at barriers of strengths in the range T a (t) < s < r a (t) + f a (t)dt are dislodged 
and move forward until they reach barriers of strength s > r a (t) + f a (t)dt. This motion of dislocations 
results in an net increase in the plastic deformation. The dislocation jumps between obstacles are assumed 
to be instantaneous. This idealization is justified when the duration of the flights is much smaller than the 
characteristic time of variation of the loads. 



As noted by Ortiz and Popov Q59Q, the information needed to describe the dislocation motion is fully 
contained in the probability density function f a (s,t), which represents the fraction of dislocation length 
facing obstacle pairs of strength s at time t. The function f a (s,t) evolves in time due to the process of 
redistribution of dislocation line described above. Initially, though, the dislocations may be assumed to be 
randomly distributed over their slip plane, and f a (s,0) = f a (s,0). At later times, f a (s,t) must vanish 
identically for < s < r a (t) in the rate independent limit. 

A kinetic equation for the evolution of f a (s, t) was derived by Ortiz and Popov [^] using standard tools 
of non-equilibrium statistical mechanics. Cuitino and Ortiz [Q] were able to obtain an analytical solution to 
this kinetic equation for the case of monotonic loading and an arbitrary time variation of the density of point 
obstacles. Remarkably, the solution takes the simple closed form 

f a (s,t) = C^- H(s-T a (t)) (34) 

Solution (Q) implies that, under the conditions of the analysis, the probability density f a (s,t) remains 
proportional to f a (s, t) over the current admissible range [r a (t), s max ]- 

Let p a (t) denote the current dislocation length per unit volume for the slip system a. The dislocation 
density released during an increment of the resolved shear stress from r a (t) to r a (t) + f a (t)dt gives rise to 



an incremental plastic strain [p8, 11] 



drf(t) = bp a (t)f (r a (t),t)) f a (t)dtN(t)l(t) (35) 

where l(t) is the average distance between obstacles, and N(t) is the average number of jumps the disloca- 
tion segments make before attaining stable positions. 

To a good approximation, l(t) can be identified with the average distance between point obstacles. As- 
suming that the obstacles are randomly distributed over the slip plane with density n a (t), a straightforward 
derivation from (17) gives ^ 



m=m)= ^m ,36) 

Next we compute the average number of jumps N(t). Evidently, the probability that an unstable segment 
becomes arrested after the first jump is equal to the probability that the first barrier encountered is of a 
strength s > T a (t). This probability is 1 — P a (T a (t),t). The probability that the segment goes beyond 
the first barrier is P a {r a {t),t). Likewise, the probability that a segment gets arrested at the second barrier 
encountered is P a (T a (t),t) [1 - P a (r a (t), t)], and the probability that it goes beyond is P a2 (r a (t),t), 
and so on. Hence, the average number of jumps between barriers taken by an unstable segment is 



N(t) = [1-P a (r a (t),t)] + 2P a (r a (t),t) [1 - P a (r a (t),t)] 
+ 3P a2 (T a (t),t) [1 - P a (r a (t),t)] + . . 



1 (37) 



1 - P a (T a (t),t) 
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Interestingly, if r Q = 0, any moving dislocation segment is sure to be arrested at the first obstacle it en- 
counters, and N = 1, in agreement with (37). Likewise, if r a {t) > s max » then the segment never reaches a 
stable barrier and N —> do, as predicted by (|37|). 



(38) 



Substituting @ and @ with s = r a (t) into @ we finally obtain 

f a (r a (t),t) 



' W ' c W [1-P«(r«(t),t)] 2 
where the characteristic plastic strain 7"(t) is defined as 

b P a (t) 



f a (t) 



7c Q (0 



Equation (38 ) defines a relation of the form 



i a {t) 



_ r a (t) 



where 



h a (t) 



h"{t) 

1 [l-P-(r"(t),t)] 2 



7c «(t) f<*(r<*(t),t) 



(39) 



(40) 



(41) 



is the hardening modulus of the slip system a. By way of example, in the particular case of obstacles of 
uniform strength, s a ^ = s max , the hardening modulus takes the form: 



h a (t) = Kit) 



T a {t)-T%\ 2 1 



7T 



cxp 



r a (t)-rf 



where 



Pq (^maxi t) 



r/edge T a (t) 

and kg®- cW 



Smax)) (42) 



(43) 



6 — 7 -(t) 

are a characteristic shear stress and plastic modulus, respectively. Equation ( |42| ) predicts an initial infinite 
hardening modulus at r = rp. The hardening modulus subsequently decreases monotonically to zero as r 
approaches s max . 

As shown in [|7p, the above relations can be extended simply to account for elastic unloading, provided 
that T a (t) does not change sign at any time during the loading history. This extension consists in defining 
the current flow stress g a (t) of slip system a as the maximum previously attained value of r a {t). The 
hardening relations are then rewritten in the form 

9 a (t) 



and 



h a {t) 



1 [1-P a (g a (t),t)} 2 
7 c a (*) f a (g a (t),t) 



(44) 



(45) 



It should be noted that, since g a (t) = whenever either r a (t) < g a (t) or f a (t) < 0, relation (44) implicitly 
accounts for elastic unloading. To account for the initial distribution of dislocations and obstacles, we take 

g a (t) = go- 



n 



5 Dislocation intersections 



In this section we proceed to estimate the obstacle strengths which enter in the preceding analysis of forest 
hardening. The interaction between primary and secondary dislocations may result in a variety of reaction 
products, including jogs and junctions [52, 53, 49|, pL 46, |], 55, 51]. Experimental estimates of junction 
strengths have been given by Franciosi and Zaoui $56(\ for the twelve slip systems belonging to the family of 
{111} planes and [110] directions in fee crystals, and by Franciosi [|7|] for the twenty-four systems of types 
{211} [111] and {110} [111] in bec crystals. The strength of some of these interactions has recently been 
computed using atomistic and continuum models [|9[ 5(J |51J]. Tang et al. have numerically estimated the 
average strength of dislocation junctions for Nb and Ta crystals po]]. 

For purposes of the present theory, we specifically concern ourselves with short-range interactions be- 
tween dislocations which can be idealized as point defects. For simplicity, we consider the case in which 
each intersecting dislocation acquires a jog. The energy of a pair of crossing dislocations is schematically 
shown in Fig. 4 as a function of some convenient reaction coordinate, such as the distance between the dislo- 
cations. The interaction may be repulsive, resulting in an energy barrier, or attractive, resulting in a binding 
energy, Fig. 4. In the spirit of an equilibrium theory, here we consider only the final reaction product, cor- 
responding to a pair of jogged dislocations at infinite distance from each other, and neglect the intermediate 
states along the reaction path. In addition, we deduce the strength of the obstacles directly from the energy 
supply required to attain the final state, i. e. the jog-formation energy. Despite the sweeping nature of these 
assumptions, the predicted saturation strengths in multiple slip are in good agreement with experiment (cf 
Section [7]), which lends some empirical support to the theory. 

We estimate the jog formation energy as follows. Based on energy and mobility considerations already 
discussed, we may expect the preponderance of forest dislocations to be of screw character, and the mobile 
dislocation segments to be predominantly of edge character. We therefore restrict our analysis to intersec- 
tions between screw and edge segments. The geometry of the crossing process is schematically shown in 
Fig. 5. Each dislocation acquires a jog equal to the Burgers vector of the remaining dislocation. The energy 
expended in the formation of the jogs may be estimated as 



pjog 



\bP\ \u aP -U aa cos8 al3 



+ \b° 



jjf3a _ jj/3/3 



COS I 



(46) 



Here, a designates the slip system of the moving edge segment, (3 the slip system of the forest screw 
dislocation, U aa = U edge is the energy per unit length of an edge segment in the slip system a, U^" = 
[7 screw is the energy per unit length of a screw segment in slip system 0, U a/3 is the energy per unit length 
of a segment in slip system a aligned with b 13 , U@ a is the energy per unit length of a segment in slip system 
[3 aligned with b a , 9 a ^ is the angle between the edge direction in system a and b^, and 6@ a is the angle 
between the screw direction in system (3 and b a . 

We additionally assume that the screw orientation defines a cusp in the dependence of the dislocation 
line energy on segment orientation, and that small deviations of a segment from a pure screw character raise 
its energy to a level comparable to U edge . This gives the energy estimate 



V screw if b a = b lS 
[/edge otherwise 



(47) 



Inserting these energies into ( |46| ) gives 



K7 



^jyscrew 
y^j screw 



[l - r cos 8 a P] 
[2r - cos(6^ a ) 



r cos 



if b a = b 13 
9 a ^l otherwise 



(48) 
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A2 A2'A3 A3'A6 A6'B2 B2"B4 B4'B5 B5'C1 C1'C3 C3"C5 C5"D1 D1"D4 D4"D6 D6' 



A2 

A2' 

A3 

A3' 

A6 

A6' 

B2 

B2" 

B4 

B4' 

B5 

B5' 

CI 

CI' 

C3 

C3" 

C5 

C5" 

Dl 

Dl" 

D4 

D4" 

D6 

D6" 



— 1.0 1.0 1.0 1.0 1.0 1.5 1.5 
1.0 — 1.0 1.0 1.0 1.0 3.23.2 
1.0 1.0 — 1.0 1.0 1.0 2.4 2.4 
1.0 1.0 1.0 — 1.0 1.0 1.8 1.8 
1.0 1.0 1.0 1.0 — 1.0 2.42.4 
1.0 1.0 1.0 1.0 1.0 — 1.8 1.8 
1.5 1.5 1.5 1.5 1.5 1.5 — 1.0 
3.2 3.2 3.2 3.2 3.2 3.2 1.0 — 
2.4 2.4 2.4 2.4 2.4 2.4 1.01.0 
1.8 1.8 1.8 1.8 1.8 1.8 1.01.0 

2.4 2.4 2.4 2.4 2.4 2.4 1.01.0 
1.8 1.8 1.8 1.8 1.8 1.8 1.01.0 
1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 
1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 

1.5 1.5 1.5 1.5 1.5 1.5 2.42.4 
3.2 3.2 3.2 3.2 3.2 3.2 1.8 1.8 
2.4 2.4 2.4 2.4 2.4 2.4 1.5 1.5 
1.8 1.8 1.8 1.8 1.8 1.8 3.23.2 
1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 
1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 

2.4 2.4 2.4 2.4 2.4 2.4 1.5 1.5 
1.8 1.8 1.8 1.8 1.8 1.8 3.23.2 

1.5 1.5 1.5 1.5 1.5 1.5 2.42.4 
3.2 3.2 3.2 3.2 3.2 3.2 1.8 1.8 



1.51.5 1.51.5 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 
3.23.2 3.23.2 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 
2.42.4 2.4 2.4 1.5 1.5 1.5 1.5 1.5 1.5 2.4 2.4 2.4 2.4 2.4 2.4 
1.8 1.8 1.8 1.8 3.2 3.2 3.2 3.2 3.2 3.2 1.8 1.8 1.8 1.8 1.8 1.8 
2.42.4 2.4 2.4 2.4 2.4 2.4 2.4 2.42.4 1.5 1.5 1.5 1.5 1.5 1.5 
1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 3.2 3.2 3.2 3.2 3.2 3.2 
1.01.0 1.01.0 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 
1.01.0 1.01.0 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 
— 1.0 1.01.0 2.4 2.4 2.4 2.4 2.4 2.4 1.5 1.5 1.5 1.5 1.5 1.5 
1.0 — 1.0 1.0 1.8 1.8 1.8 1.8 1.8 1.8 3.2 3.2 3.2 3.2 3.2 3.2 
1.01.0 — 1.0 1.5 1.5 1.5 1.5 1.5 1.5 2.4 2.4 2.4 2.4 2.4 2.4 
1.01.0 1.0 — 3.2 3.2 3.2 3.2 3.2 3.2 1.8 1.8 1.8 1.8 1.8 1.8 
1.8 1.8 1.8 1.8 — 1.0 1.0 1.0 1.0 1.0 3.2 3.2 3.2 3.2 3.2 3.2 
1.8 1.8 1.8 1.8 1.0 — 1.0 1.0 1.0 1.0 3.2 3.2 3.2 3.2 3.2 3.2 
2.42.4 2.4 2.4 1.0 1.0 — 1.0 1.0 1.0 2.4 2.4 2.4 2.4 2.4 2.4 
1.8 1.8 1.8 1.8 1.0 1.0 1.0 — 1.0 1.0 1.8 1.8 1.8 1.8 1.8 1.8 
1.5 1.5 1.5 1.5 1.0 1.0 1.0 1.0 — 1.0 2.4 2.4 2.4 2.4 2.4 2.4 
3.23.2 3.23.2 1.0 1.0 1.0 1.0 1.0 — 1.8 1.8 1.8 1.8 1.8 1.8 
1.8 1.8 1.8 1.8 3.2 3.2 3.2 3.2 3.2 3.2 — 1.0 1.0 1.0 1.0 1.0 
1.8 1.8 1.8 1.8 3.2 3.2 3.2 3.2 3.2 3.2 1.0 — 1.0 1.0 1.0 1.0 
1.5 1.5 1.5 1.5 2.4 2.4 2.4 2.4 2.4 2.4 1.0 1.0 — 1.0 1.0 1.0 
3.23.2 3.23.2 1.8 1.8 1.8 1.8 1.8 1.8 1.0 1.0 1.0 — 1.0 1.0 
2.42.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 2.4 1.0 1.0 1.0 1.0 — 1.0 
1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.8 1.0 1.0 1.0 1.0 1.0 — 



Table 2: Normalized jog-formation energies resulting from crossings of bcc dislocations. 



where r = [/edge /[/screw j s t ^ e rat j Q Q f screw to ec ig e dislocation line energies. For Ta, recent atomistic 
calculations [|^] give a value of r = 1.77. The resulting jog formation energies for the complete collection 
of pairs of {211} and {110} dislocations are tabulated in Table|2[ 



The net effect of jog formation on hardening may be ascertained as follows. Consider the case in which 
a primary system a contains forest obstacles of a single species, corresponding to secondary dislocations 
in the slip system j3. Let n a " be the density of forest obstacles per unit area of the primary plane a. As a 
primary dislocation sweeps through an area A, the energy expended in forming jogs with forest dislocations 
of the /3-type is E } °^n al3 A. On the other hand, the potential energy released as a result of the motion of the 
dislocation follows from the Peach-Koehler formula as br a A, where r a is the resolved shear stress acting on 
the primary system a. Hence, the forest obstacles of type (3 become 'transparent' to the motion of primary 
dislocations when br a A > E 3 °§n a ^A, or 

r a > \E^n^ (49) 

Since the jog energies scale with the elastic moduli, they may be expected to reduce to zero at the melting 
temperature, which to some extent accounts for the observed thermal softening. In addition, the jog energies 
are small enough for thermal activation to be operative at the level of individual obstacles. A derivation 
entirely analogous to that leading to Eq. (Q) yields, in this case, 

s a0 1 . , (T (3 

— ^ = i — asmh — e 
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where 



af = - af3 (51) 

fta £junct 



and 

7« = 2p a bl a u D (52) 

The lengths T a and l) nnct describe the geometry of the junction, as illustrated in Fig. 6, and are of the order 
of a few Burgers vectors. 

It follows from these developments that, at low temperatures, the energies collected in Table ||| provide 
a measure of the corresponding obstacle strengths. We may also recall that the saturation flow stress of a 
slip system is directly related to the strengths of the obstacles. Consequently, information about obstacle 
strengths may be inferred from saturation flow stresses in crystals deformed in double slip. Measurements 



of this type have been performed by Franciosi for a-Fe J27|]. The general trends exhibited by Table g are 
consistent with Franciosi's data, as well as with the recent analytical model of Lee et al. [|67|]. In addition, the 
saturation flow stresses resulting from the model developed above are in good quantitative agreement with 
those measured experimentally, cf Section ^, which lends empirical support to jog-formation as a plausible 
mechanism underlying short-range obstacle strength and saturation in bcc crystals. 



6 Dislocation evolution 

The density of forest obstacles depends directly on the dislocation densities in all slip systems of the crystal. 
Therefore, in order to close the model we require a equation of evolution for the dislocation densities. 
Processes resulting in changes in dislocation density include production by fixed sources, such as Frank- 



Read sources, breeding by double cross slip and pair annihilation (see Q68Q for a review; see also Q69| , |70| , [71 



72, [73|, [74]]). Although the operation of fixed Frank-Read sources is quickly eclipsed by production due to 
cross slip at finite temperatures, it is an important mechanisms at low temperatures. The double cross-slip, 
fixed Frank-Read sources and pair annihilation mechanisms are next considered in turn. 

6.1 Breeding by double cross-slip 

The importance of breeding by double cross-slip as a dislocation generation mechanism in crystals was 
emphasized by Johnston and Gilman [69], 7(J. In this mechanism, the screw section of a moving dislocation 



migrates to a parallel plane by double cross-slip, thus creating a pair of sessile segments which pin down the 
dislocation and act in a manner similar to a Frank-Read source. 

The rate at which dislocation length is produced by this mechanism can be estimated as follows. Let 
Nq S denote the number of dynamic sources per unit volume induced by double cross-slip on slip system a. 
Additionally, let L sat denote the dislocation length emitted by each source prior to saturation. The rate of 
dislocation production per unit volume is, therefore: 

pg s = L^Ng s (53) 



Assuming that cross-slip is thermally activated (see, e. g., J751, |7q] for recent calculations of pathways and 
energy barriers in copper), the rate Nq S at which dynamic sources are generated may in turn be computed 
as follows. Let L be the mean-free path between cross-slip events, L cross the length of the screw segment 
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effecting the double cross glide, and E cross the energy barrier for cross slip. The frequency of cross glide at- 
tempts per unit volume is, therefore, p a v a / {LL CTOSS ), where v a is the mean dislocation velocity. According 
to transition-state theory, of these attempts a fraction e _ ^ £;cross is successful. This gives: 

Ng s = Pl_ e -PB— (54) 

7* rcross v ' 

Indeed, theoretical J77| ] and experimental [69, 7(J investigations suggest that the breeding rate due to cross- 
glide is proportional to the mean dislocation speed. Using Orowan's formula, ( p4| ) may be rewritten as 

Arg = JF—e-PE?*- (55) 
bLL cross 



Inserting this relation into (53) gives 



rsat ■ a 

•a ^ 7 „-/3E CIOSS , <f -s 
Prs = T"TT e (JO) 

bLL CTOSS 



Following Essmann and Rapp [72], we shall additionally assume that the mean free path is inversely pro- 



portional to the dislocation density, which gives: 

bph = XcsVp^i a (57) 

where we have written 

r sat 

Acs = i_ e -^ cross (58) 

^ J rcross v ' 

An improvement on this model is to account for the resolved shear stress acting on the cross-slip plane in 
the computation of the activation energy, resulting in a so-called non-Schmid effect, but this possibility will 
not be pursued here in the interest of simplicity. 

The double cross-slip mechanism is geometrically similar to the double kink formation alluded to earlier. 
By virtue of this similarity, we may expect that L cross m L and E cvoss m E kink to a first approximation. 

6.2 Multiplication by fixed Frank- Read sources 

An infinite amount of dislocation multiplication can be sustained by the fixed Frank-Read sources as long 
as newly formed loops can expand and move away from the pinning points. The precise details of this 



mechanism, which was independently proposed by Frank and Read[|78|] in the early 50, are well documented 



and can be consulted elsewhere[39]. In here, we limit our analysis to obtained an estimate for the dislocation 
multiplication based on this mechanism. Considering that the obstacles in the slip plane a serve as pinning 
points to operate the Frank-Read sources, the increase of the dislocation population is given by 

Pfr = LfrN^r (59) 
where 7Vp R is the number of active sources per unit volume and Lp R is the length of the emitted segment 



by the source and n a is the obstacle density given by Eq. (|25j). The effective rate of activation of Frank- 
Read sources in a given system a can be estimated by the number of intersections of the gliding loops with 
obstacles in that system weighted by the efficiency of the source. In order to account for the possibility of 
the newly formed loops to glide away from the obstacles, we introduce the factor \J p a /n a which is the ratio 
between mean distance between the obstacles and the mean dislocation spacing. Then, 

N^ R = J Pl n a v a p a . (60) 
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Introducing Orowan's formula and assuming that the length of dislocation segments is proportional to the 
mean distance between obstacles, the multiplication rate due to fixed Frank-Read sources can be written as 



bp a 



FR 



Afr\/p^7 c 



(61) 



Unlike breeding by cross-glide, dislocation multiplication by fixed Frank-Read sources is not thermally 
activated process, and thus, it remains operative even at low temperatures. At finite temperatures however 
both mechanisms contribute to the total dislocation multiplication rate. 

It should pointed out that both Eqs. (57) and (^Tj) predict that p a grows as (j 01 ) 2 - This rate of growth is 
indeed observed in many crystals during the intermediate stages of dislocation multiplication. 



6.3 Attrition by pair annihilation 



Dislocation densities are often observed to attain a saturation density at sufficiently large strains Q79Q . This 
saturation stage arises as a result of the competition between dislocation multiplication mechanisms such as 
double cross-slip and pair annihilation. For instance, spontaneous annihilation is observed in metals such as 



copper for screw dipole heights less than 1 nm [|80p. Pair annihilation is mainly the result of the cross-slip 
of screw segments of opposite sign 



1 ,35,621 



Huang et al. have studied the dynamic stability of short-range linear-elastic interactions between 
two dislocations of parallel line vectors which glide on two parallel slip planes in bcc crystals. Here we 
develop a similar but somewhat simpler linear-elastic model of dislocation pair annihilation. Thus, our goal 
is to estimate the frequency with which two parallel screw segments moving on parallel planes will acquire 
converging trajectories leading to their mutual annihilation. We note that, for simplicity, our analysis is 
restricted to annihilation between pairs of dislocations belonging to the same slip system. 

Consider a screw segment moving under the action of an applied shear stress r. The segment follows a 
path which brings it in close proximity to a second immobile screw segment, Fig. 7. Let (x, y) be coordinates 
centered at the fixed dislocation such that x points in the direction of motion of the incoming dislocation 
and y points in the direction of the cross-slip plane, Fig. 7. The interaction force per unit length exerted on 
the moving segment along the cross slip plane is [^9|] : 



fv 



-K 



y 



where K is the pre-logarithmic factor for a screw segment, and 

r = \J x 2 + y 2 + 2xy cos 



o 



(62) 



(63) 



where a is the angle subtended by the x and y directions, Fig. 7. For bcc crystals, a = 2ir/3. For isotropic 
crystals, K = ph 2 /4ir. In view of Eq. (11), the equation of motion for the incoming screw dislocation is: 



y 



a 771 cross 

-2l P u D e~ pE sinh 



(3E C 



, (Ky/br 2 

TO 



for y > kq > 0, where the characteristic distance 



K 



K 
br 



(64) 



(65) 



is determined by the condition f y = 6ro- For y < kq dislocation motion is no longer a thermally ac- 
tivated process and therefore Eq. |64], which is based on the transition-state theory, does no longer ap- 
ply. In this regime, where the interaction force exceeds the Peierls barrier, dislocation mobility increases 
sis nificantly[|47]]. 
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Depending on the value of y as x — > — oo, the moving segment bypasses or is captured by the second 
segment. The precise calculation of the trajectories of the moving dislocation requires numerical computa- 
tion. However, the essence of the annihilation mechanism can be captured by the following simple argument. 
Consider a dislocation initially at rest at position (x = 0, yo)- Subsequently, the motion of the dislocation 
for y > Ko is governed by the equation: 

y = -2l P u D e-^° ss sinh (pi^!£M\ (66) 

This equation is separable, and the time required for the annihilation process follows simply as: 

t = | 2 /p^e-^ cross sinh (^E^^M^y 1 dy (67) 

In order to facilitate this calculation, the time required to travel from y = kq to y = is neglected in the 
previous equation due to the much higher dislocation mobility in this regime and the small value of kq, of 
the order of few b for BCC crystals. Additionally, we may simply set sinh(x) x to a first approximation, 
with the result: 

2l P v D f3E CIOSS 2k q 

On the other hand, the time which the incoming dislocation spends at distances of the order of yo to the 
receiving dislocation is 

t ~ -^2/o (69) 



For annihilation to be possible the annihilation time (68) must be less than the time (59) which the incoming 



dislocation spends in the immediate vicinity of the target dislocation. This yields the condition: 

yo < k = k (A + VA 2 + 1) (70) 

where 

A = e -/3£ cr ° ss / 3£cross^cross^a (71) 

is a factor depending on the strain rate and temperature, 

7 cross = 2bpl P u D (72) 

is a reference slip-strain rate, and k may now be regarded as an effective pair annihilation distance. The 
cut-off value k c , corresponding to the maximal pair annihilation distance, is the effective screening distance 
which can be set equal to the mean distance between dislocations. 

A simple expression which interpolates between the extreme values of k is 

1 1 1 

- = - + — ( v (73) 

K c k La + va 2 + rj 

It follows that the critical pair-annihilation distance k decreases with increasing strain rate and decreas- 
ing temperature. Thus, at high strain rates the dislocation velocities are high and the probability of being 
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captured by another dislocation diminishes accordingly. Additionally, an increase in temperature increases 
the dislocation mobility and speeds up the annihilation process, which results in an attendant increase in 
annihilation rates. As will be demonstrated in Section fTj, these trends are important in order to capture the 
temperature and strain-rate dependence of the stage I-II transition in Ta. 

The rate of dislocation attrition due to pair annihilation may finally be estimated as follows. Over a 
differential of time dt, the number of annihilation events per unit volume is 



dN a 



1 



( P a ) 2 KV a dt 



2L C 



The attendant loss of dislocation length per unit volume is 



dp° 



-2L cross dN c 



(74) 



(75) 



Combining these relations and using Orowan's formula, the rate of pair annihilation may finally be expressed 
in the form 



-np a j a 



(76) 



6.4 Dislocation multiplication rate 



Combining (|57|), (pl|)) and (|76|), the total rate of change of the dislocation density follows as 

bp a = \^j a - K p a j a 

where A = Acs + ^fr- Evidently, p a = upon the attainment of the saturation density: 



(77) 



„sat 



Dividing through by p sat , (fTJ) may be recast in the form: 



£_ 

,,sat 



where 



(78) 



(79) 



,,sat 



(80) 



is a saturation slip strain. Equation ( |79| ) can be integrated, yielding: 



sat 




(81) 



The rate equation (175) expresses a competition between the dislocation multiplication and annihilation 



mechanisms. For slip strains j a <C 7 
dislocation density p a grows as (7 ) 2 



Silt 



, the multiplication term dominates and, as noted previously, the 
By contrast, when 7" 3> 7 sat the rates of multiplication and an- 
nihilation balance out and saturation sets in. After saturation is attained, the dislocation density remains 



essentially unchanged. It should be carefully noted, that, in view of <J73|), the saturation slip strain 7 sat is 
a function of temperature and strain rate. In particular, 7 sa t decreases with increasing temperature and de- 
creasing strain rate. Since the stage I-II transition strain scales roughly with 7 sat , we expect these trends to 
be exhibited by the transition strain itself, in accordance with experimental observation. 
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Table 3: Parameter set for Tantalum 



Parameter 


Value 


Units 


Ch 


266.49 - 0.021T (*) 


[GPa] 


Cl2 


156.25 - 0.006T 


[GPa] 


C44 


90.02 - 0.015T 


[GPa] 


b 


2.86 x 10" 10 


[m] 


,77^ kink 


0.70 


[eV] 


L kink /6 


13 




[/ edge //z6 2 (**) 


0.216 




edge / screw 


1.77 




l/b 


5 




L Junct /6 


20 




^cross 


0.67 


[eV] 


L CTOaa /b 


13 




K c /b 


1250 




Afr 


2.3 




Po 


10 12 


[m- 2 ] 


go 


8.0 


[MPa] 


a 


0.01 





* T in Kelvin. 



** M — 1^44 + g(Cll — C12). 



7 Comparison with experiment 

We proceed to validate the theory against the uniaxial tests on Ta single crystals of Mitchell and Spitzig 
In these tests, 99.97%-pure Ta specimens were loaded in tension along the [213] crystallographic axis, at 
various combinations of temperature and strain rate. In particular we considered temperatures ranging from 
296 K to 573 K, and strain rates ranging from 10 _1 s _1 to 10~ 5 s _1 . The numerical procedure employed 



for the integration of the constitutive equations has been described elsewhere [g8J. The constitutive update is 
fully implicit, with the active systems determined iteratively so as to minimize an incremental work function. 
All stress-strain curves are reported in terms of nominal stress and engineering strain. 

The material property set used in calculations is collected in Table The elastic moduli C\\, C\i 
and C44 were obtained by fitting to the tables of Simmons and Wang p2|]. The ratio between edge and 
screw dislocation-line energies (f/ ed g e /[/ screw ) j s taken from the atomistic calculations of Wang et al. [f35]. 
The remaining parameters have been obtained by fitting to the experimental data of Mitchell and Spitzig [g]. 
Note that we do not give any value for L sat , since it has no measurable influence in the range of temperatures 
considered here. It is expected, however, to play a role at higher temperatures. 

Figs. 8 and 9 show the predicted and measured stress-strain curves for a [213] Ta crystal over a range 
of temperatures and strain rates. It is evident from these figures that the model captures salient features of 
the behavior of Ta crystals such as: the dependence of the initial yield point on temperature and strain rate; 
the presence of a marked stage I of easy glide, specially at low temperature and high strain rates; the sharp 
onset of stage II hardening and its tendency to shift towards lower strains, and eventually disappear, as the 
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temperature increases or the strain rate decreases; the parabolic stage II hardening at low strain rates or high 
temperatures; the stage II softening at high strain rates or low temperatures; the trend towards saturation at 
high strains; and the temperature and strain-rate dependence of the saturation stress. 

The theory reveals useful insights into the mechanisms underlying these behaviors. For instance, since 
during state I the crystal deforms in single slip and the secondary dislocation densities are low, the Peierls 
resistance dominates and the temperature and strain-rate dependency of yield owes mainly to the thermally 
activated formation of kinks and crossing of forest dislocations. It is interesting to note that during this stage 
the effect of increasing (decreasing) temperature is similar to the effect of decreasing (increasing) strain rate, 
as noted by Tang et al. $6&. The onset of stage II is due to the activation of secondary systems. The rate 
at which these secondary systems harden during stage I depends on the rate of dislocation multiplication in 
the primary system. This rate is in turn sensitive to the saturation strain 7 sat , which increases with strain 
rate and decreases with temperature. As a result, the length of the stage I of hardening is predicted to 
increase with strain rate and decrease with temperature, as observed experimentally. Finally, the saturation 
stress is mainly governed by the forest hardening mechanism and, in particular, by the strength of the forest 
obstacles. This process is less thermally activated than the Peierls stress, since the corresponding energy 
barriers are comparatively higher. Consequently, the stress-strain curves tend to converge in this regime, in 
keeping with observation. 

The apparent softening observed in simulation results at the lowest temperature (296 K) and the highest 
strain rate (1CP 1 s _1 ) is actually an effect of the boundary conditions, which allow for some level of rotation 
of the specimen. Since in those cases, the material hardening is relatively low (stage I only), this geometrical 
softening dominates the apparent macroscopic behavior. In the other cases, the activation of several systems 
limits the extend of the rotations reducing the effects of the macroscopic hardening. In order to simulate 
more precisely the experimental boundary conditions, a model of the entire specimen allowing for a non- 
homogeneous deformation field should be considered. 

The effect of temperature and strain-rate on hardening is also illustrated by the evolution of slip activity 
in the primary and secondary slip systems, Dl and DA", respectively. Figs. 10 and 11 show the evolution 
of slip strains and dislocation densities as a function of temperature. It is evident from these figures that 
dislocation density saturation in the primary system occurs earlier as the temperature is increased, resulting 
in the activation of the secondary system and in the onset of stage II at lower strains. A similar effect is 
observed when the strain rate is decreased, Figs. 12 and 13. 

The effect of loading direction on the hardening rate is illustrated in Fig. 14. The figure shows the stress- 
strain curves obtained by loading the crystal in the [213], [101] and [111] directions at a strain rate of 10~ 3 
s _1 and a temperature of 373 K. Each of these loading directions results in the activation of a different set 
of slip systems. As may be observed in the figure, the higher initial yield stress in the [111] direction 
relative to the baseline [213] direction, and the initial negative hardening rate, are fairly well captured 
by the model. The experimental curve exhibits a subsequent upturn, most likely due to the activation of 
additional secondary systems, which is not captured by the model. Loading in the [101] direction results 
in the activation of a large number of systems from the outset and the rate of hardening is correspondingly 
high. The model appears to over-predict the rate of hardening. 

However, as was noted in [0], the stress-strain curve for crystals loaded in high-symmetry orientations 
are extremely sensitive to small misalignments in the loading axis, which accounts for the large experimental 
scatter characteristic of those orientations. This extreme sensitivity is due to the fact that small deviations 
from a high-symmetry loading axis, of the order of a degree or less, may result in the activation of a different 
set of slip systems, which may in turn have a large effect on the hardening rates. This pathological behavior 
is illustrated in Fig. 15, which shows the stress-strain curve obtained by randomly offsetting less than one 
degree the nominal loading axis [101]. The resulting reduction in the hardening rate is quite remarkable. 
In addition, the experimental curve falls between the bounds of the stress-strain curves for the nominal 
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and perturbed directions. This comparison exemplifies the need to take experimental scatter into account 
when assessing the fidelity of models, specially where crystals loaded along directions of high symmetry 
are concerned. 

8 Summary and conclusions 

We have developed a micromechanical model of the hardening, rate-sensitivity and thermal softening of bcc 
crystals. The model is predicated upon the consideration of an 'irreducible' set of unit processes, consisting 
of: double-kink formation and thermally activated motion of kinks; the close-range interactions between pri- 
mary and forest dislocation, leading to the formation of jogs; the percolation motion of dislocations through 
a random array of forest dislocations introducing short-range obstacles of different strengths; dislocation 
multiplication due to breeding by double cross-slip; and dislocation pair-annihilation. Each of these pro- 
cesses accounts for-and is needed for matching-salient and clearly recognizable features of the experimental 
record. In particular, on the basis of detailed comparisons with the experimental data of Mitchell and Spitzig 
the model is found to capture: the dependence of the initial yield point on temperature and strain rate; 
the presence of a marked stage I of easy glide, specially at low temperature and high strain rates; the sharp 
onset of stage II hardening and its tendency to shift towards lower strains as the temperature increases or the 
strain rate decreases; the initial parabolic hardening followed by saturation within the stage II of hardening; 
the temperature and strain-rate dependence of the saturation stress; and the orientation dependence of the 
hardening rates. 

The choice of analysis tools which we have brought to bear on the unit processes of interest, e. g., 
transition-state theory, stochastic modeling, and simple linear-elastic models of defects and their interac- 
tions, is to a large extent conditioned by our desire to derive closed-form analytical expressions for all 
constitutive relations. As noted throughout the paper, many of the mechanisms under consideration are 
amenable to a more complete analysis by recourse to atomistic or continuum methods. However, at this 
stage of development, direct simulation methods, be it atomistic or continuum based, tend to produce un- 
manageable quantities of numerical data and rarely result in analytical descriptions of effective behavior. 
The daunting task of post-processing these data sets and uncovering patterns and laws within them which 
can be given analytical expression is as yet a largely unfulfilled goal of multiscale modeling. 

This larger picture notwithstanding, one concrete and workable link between micromechanical models 
and first-principles calculations concerns the calculation of material constants. A partial list relevant to the 
present model includes: energy barriers and attempt frequencies for double-kink formation, kink migra- 
tion, dislocation unpinning, cross-slip, and pair annihilation; dislocation-line and jog energies; and junction 
strengths. Other properties which have yielded to direct calculation include the volumetric equation of 
state (EoS), the pressure dependence of yield, and the pressure and temperature dependence of elastic mod- 
uli. References to recent work concerned with the calculation of these material properties have been given 
throughout the paper. As noted earlier, these results provide a suitable basis for future extensions of the 
present model to higher temperatures, pressures and strain-rates. 
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